Response of a Hodgkin-Huxley neuron to a high-frequency input 
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We study the response of a Hodgkin-Huxley neuron stimulated by a periodic sequence of conduc- 
tance pulses arriving through the synapse in the high frequency regime. In addition to the usual 
excitation threshold there is a smooth crossover from the firing to the silent regime for increas- 
ing pulse amplitude Qsyn- The amplitude of the voltage spikes decreases approximately linearly 

with Qsyn. 

In some regions of parameter space the response is irregular, probably chaotic. In the chaotic 
regime between the mode- locked regions 3:1 and 2:1 near the lower excitation threshold the output 
interspike interval histogram (ISIH) undergoes a sharp transition. If the driving period is below 
the critical value, Ti < T*, the output histogram contains only odd multiples of Ti. For Ti > T* 
even multiples of Ti also appear in the histogram, starting from the largest values. Near T* the 
ISIH scales logarithmically on both sides of the transition. The coefficient of variation of ISIH has 
a cusp singularity at T*. The average response period has a maximum slightly above T*. Near the 
excitation threshold in the chaotic regime the average firing rate rises sublinearly from frequencies 
of order 1 Hz. 



I. INTRODUCTION 

Biological neurons transmit information in the form of 
sharp spikes of potential difference across the lipid bi- 
layer forming the wall of the nerve cell. This feature of 
the cell's reaction to input signals is remarkably consis- 
tent in different organisms and different types of neurons. 
The action potential spikes are assumed to be the prin- 
cipal carrier of information. The early view that infor- 
mation is transmitted via rate coding has evolved. It is 
now recognized that also the spike time coding is used 
in neural systems^'^. While the precise coding recipe is 
unknown it is clear that the knowledge of the response 
of various types of neurons to different stimuli is funda- 
mental to formulating the theory of information transfer 
in the neural system. 

Our understanding of conductance-based models of 
neurons is largely based on the Hodgkin-Huxley (HH) 
model originally formulated to describe the dynamics of 
the membrane potential of the squid giant axon^. The de- 
tailed voltage-clamp measurements of the voltage-gated 
potassium and sodium ion currents led to revisions of the 
HH model. The modifications required to achieve bet- 
ter agreement with experiments were reviewed by Clay^. 
Studies of single neurons and neuronal networks often 
employ simplified models, such as integrate- and- fire and 
FitzHugh-Nagumo (FHN) models^^^. It is believed that 
the two-dimensional flow models such as FHN reproduce 
qualitatively the behavior of the HH model. However 
these simplifications are not always justifiable^^. In an 
interesting analysis of chaos in the HH model Gucken- 
heimer and Oliva^ point out that even the concept of a 
firing threshold may be more subtle that just a smooth 
hypersurface dividing subthreshold and suprathreshold 
membrane potentials. 

Over the years many studies of HH equations were 
carried out, including stochastic variations of various 
quantitiesiSrJ^. An important question is to what extent 
the qualitative properties of neuron response depend on 



the functional form of the input signal. One frequently 
used form of input is constant plus a sinusoidal term. 
However the physiological signals are more pulse-like. In 
a strongly nonlinear system this may lead to substantial 
differences in the output. 

In the sinusoidally driven HH model the excitation 
threshold rises sharply at large frequencies. The phase 
diagram in the frequency-current amplitude plane con- 
sists of three phase locked regions with integer ratio of 
the output period to the input period, To/Tj, 1:1, 2:1, 
and 3:1. There are also areas of fractional locking and 
bistable or chaotic response around these phase-locked 
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It was pointed out that the edges of mode-locked 
plateaus have analogies to phase transitions in the equi- 
librium statistical mechanics. Two forms of scaling of the 
average deviation from perfect mode-locking were found 
near the edges of plateaus with constant p/q, where p 
and q are integers, indicating number of input spikes per 
number of output action potentials—. The scaling has 
either exponent 1/2 or is logarithmic. In this paper we 
will show that scaling is more common and appears also 
near the multimodal transition points. 

Here we assume the a form of postsynaptic current, 
Isyn ^ texp(— t/r), where t is time from the onset of 
the input spike and r is the time scale of the synaptic 
action. This form is close to experimental observation 
although it does not take into account a more complex 
dynamics of the ion channel kinetics, usually described 
in the Markovian scheme. 

The general form of the phase diagram of the Hodgkin- 
Huxley model with this input was studied initially in 
Ref^ii. However many important questions are still to 
be answered. One of them is the behavior of the system 
in the high-frequency limit. In the following we present 
the model and show the main features of high-frequency 
response. 



II. THE MODEL 

The Ho dgkin- Huxley neuron subject to periodic con- 
ductance pulses is defined by the following set of 
equations^ 



CdV/dt = -gNam^h{V - VNa) - 9Kn\V - Vk) 

-9l{V - Vl) + lext + Isyn, 



dm/dt = -{am + bm)m + a^, 



dh/dt = -{ah + bh)h + a/,, 



dn/dt = —{an + bn)n + a^, 
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1: Sample voltage trace for a constant input current 
= lOAtA/cm^ 



0.1(y + 40)/[l-e-(^+4°)/i°], (5) 



= 4e-(^+65)/i8 



ah = 0.07e-(^+65)/20, 



6;. = l/[l + e-(™)/i°], 



a„ = 0.01(V + 55)/[l-e 



-(y+55)/10i 



bn = 0.125e-(^+65)/80. 
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In equations (|5])- (p!Q]) the voltage is expressed in mV 
and the rate constants a and P are given in ms~^. 
The reversal potentials of sodium, potasium and leak- 
age channels are VNa = 50mV, Vr = — 77mV, and 
Vl = — 54.5mV, respectively. The corresponding max- 
imum conductances are gNa = 50mV, gx = 36mS/cm^, 
and gL = 0.3mS/cm^. The capacity of the membrane is 
C = l/iF/cm^3. 

The synaptic current Isyn is given by the following 
equation, 



Isyn{t) = gsyn^C^{t - tm)(K - Vsyn) , (H) 
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FIG. 2: The ratio of the average output spiking rate to the 
input rate, k = To/Ti. Mode-locked regions with /c = 1,2,3 
and /c = 4, 5, 6, 7, 8 are shown in black and grey respectively. 
Voltage peaks were counted as spikes when V exceeded 0. For 
high values of Qsyn the neuron does not respond. 



where tin denotes the start of the n*^ pulse, gsyn is the 
conductivity of the synapse, Va = 30m V is the maximum 



potential in the postsynaptic area and Vsyn = —^OmV is 
the reversal potential of the synapse. The period of the 
synaptic drive is Ti = tin+i — Un- The external current 
lext is set to 0, except for a sample run shown in Fig. [TJ 
The time-dependence is given by the function 



a{t) = (t/r)e-^/")6(t). 



(12) 



where r is time scale characterizing the dynamics of the 
synaptic action and Q{t) is the Heaviside step function. 
We study the dependence of the output interspike sepa- 
ration To on Ti and gsyn- 
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Equations (p!])- (p!Q|) were integrated with the fourth or- 
der Runge-Kutta scheme. The time step was 0.01 ms. 
For each parameter set the simulation was run for 30 
seconds. Results of the initial three seconds of each data 
set were discarded to avoid transient behavior. In the 
chaotic regime the data were obtained from five runs for 
each value of the horizontal coordinate. 



III. RESULTS 

The average output spiking rate in the form of a color 
map as a function of the input period Ti and max- 
imum synaptic conductivity ggyn is presented in Fig. 
[21 The mode-locked regions are shown as areas of uni- 
form color. For small Ti the total incoming current is 
approximately constant with a small modulation, and 
the excitation threshold rises linearly with increasing T^, 
gsyn — 0.04Ti mS/(ms cm ). For gsyn exceeding approx- 
imately OATi mS/(ms cm ) the spiking action does not 
occur. We can see from Fig. [2] that this behavior sets in 
below Ti 6ms . 

The obtained phase diagram is qualitatively different 
from a response to a sinusoidal input, where the excita- 
tion threshold diverges as l/T^, for Ti 0. In general 
we may expect that the constraint of charge balancing, 
j^t+T^ _ ^^^^ have a significant impact at high fre- 
quencies. For intermediate values of the input period, 
5ms < Ti < 13ms, the topology of the phase diagram 
resembles results obtained with sinusoidal input, see e.g. 
Fig. 2 of Refi^. 
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FIG. 4: For high synaptic conductivities the distinction be- 
tween action potential and the background oscillations loses 
its meaning. This sample was obtained for Ti = 4.5ms and 
gsyn — 2.35mS/cm^. 



Fig. [3] shows dependence of minima and maxima of 
V on gsyn for three input frequencies. The amplitude of 



response decreases linearly with increasing gsyn- There is 
no well-defined spiking threshold. There are intervals of 
parameter values for which the response is highly irregu- 
lar and the values of maxima and minima of V vary sig- 
nificantly. 
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FIG. 5: (a) The spectrum of interspike separations of the 
output signal as a function of the input period Ti for Qsyn = 
0.4mS/cm^, (b) Detailed view of the chaotic region between 
Ti — 5ms and 6ms. Each ISI cluster belongs to different /c, 
where k — 2,3,4,5,.... The distinction between k — 2 and 
A; = 3 is blurred. 
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FIG. 6: Scaling of the excitation edge of (a) odd-only mul- 
tiples of the input Tj, and (b) all integer multiples, in the 
chaotic region between k — 2 and k — ?>. For g — 0.2mS/cm^, 
the transition occurs at T = 6.54175ms. 



A sample time-dependence of the membrane potential 
is shown in Fig. HJ The maxima of V span almost the 
entire range between — 60mV and OmV. There is no clear 
separation of spikes from the rest of the signal. 

Chaotic behavior in the parameter space between the 
3:1 and 2:1 mode- locked regions leads to multimodal re- 
sponse. The interspike separation for Qsyn = 0.4mS/cm^ 
is shown in Fig. [5l For Ti between 5.5ms and 6ms all 
integer multiples of input Ti with the exception of the 
lowest one appear in the output ISIH. 

It is interesting to note that ISI histograms (ISIH) from 
some older experiments on nerve fibers of monkeys and 
single neurons in the primary visual cortex of a catii^ 
show some similarity to Fig. [5l Experimental histograms 



are sequences of diminishing peaks occuring at integer 
multiples of the input interspike separation. In Fig. [5] 
the lowest element of the sequence is missing due to the 
refractoriness of the neuron. Similar form of ISIH was ob- 
tained in a theoretical study of a bistable system stimu- 
lated by periodic function with additive Gaussian noise^ 
where the presence of noise was essential. However the 
multimodal histogram was also obtained in a simulation 
of a deterministic modification of the HH model^^. 

The HH model studied here does not contain stochastic 
terms. The multimodal response in Fig. [5] is a result of 
a deterministic nonlinearity. Thus noise is not the only 
ingredient enabling the reproduction of the multimodal 
experimental ISIH. It is possible to identify the source of 
multimodality by studying ISIH in more detail. 
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FIG. 7: The multimodal transition at ^ = 0.17mS/cm^. 
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FIG. 8: Coefficient of variation for g = 0.17mS/cm^. The 
variability near Ti = 6.6ms is due to the proximity to the 
firing threshold. 
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FIG. 9: The ratio k = To/Ti for g = 0.17mS/cm^ The 
maximum of k is shifted approximately 0.2 ms to the right 
relative to maximum of CV (see Fig. [8]). 
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FIG. 10: The location of the multimodal transition (filled 
squares) on the response diagram. 



Close to the excitation threshold, at gsyn — 
0.2mS/cm^, there exists a transition from the odd-only 
ISIH to ISIH with all integer multiples of Tj, see Fig. [6l 
Near the transition the edges of high-/c clusters scale loga- 
rithmically. The scaling holds both along the Ti axis and 
along the gsyn axis. It can be viewed as a competition 
between the odd and the even multiples of the driving 
period. 

A clear indication of this "spectral" transition is the 
singular behavior of the coefficient of variation, see Fig. 
[HI At the transition CV is of order 1. and k is signif- 
icantly larger than 3. The maximum k occurs approx- 



imately 0.2 ms above the singularity of CV. One may 
also think of this shift as a result of relaxation from the 
constraint of odd-only modes below T*. At T* the high- 
est even modes become available and this leads to the 
increase of k. 

If such transition were found experimentally it would 
be a clear sign of the deterministic nonlinear dynam- 
ics. In the presence of noise this sharp feature would 
be smeared and would vanish if noise dominates the dy- 
namics of the system. 
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FIG. 11: Frequency vs. synaptic conductivity for Ti — 
7.45ms. Each data point is averaged over 15 runs for 60 s 
with different initial conditions. The initial 6 s from each run 
were discarded. 



IV. CONCLUSIONS 

For high synaptic drives in the high frequency regime 
distinguishing the action potential from the background 
activity becomes problematic. In this limit the neuron is 
very sensitive to small changes of the functional form of 
the signal. For periodic drive with small time constantr 
and Ti below 6ms the width of the spiking regime along 
the Qsyn axis scales linearly with Ti. The quality of the 
neuron's response deteriorates linearly with increasing 
Qsyn- This is in contrast to findings for a sinusoidal sig- 
nal, and more generally for a class of signals satisfying the 
constraint of charge balancing, where the spiking action 
remains well defined in the high-frequency limit. 

A mechanism of suppression of the neuron's activity 
might help explain self-regulating behavior of neocorti- 
cal networks. Various mechanisms of homeostatic action 
for neural microcircuits were proposed. It would be use- 
ful to investigate whether more realistic extensions of the 
Hodgkin-Huxley model also exhibit self-regulation in re- 
sponse to high-frequency inputs. The network of such 
neurons would have a " safety switch" built in at the level 
of individual cells. For Ti between 4 and 6 ms the up- 
per critical synaptic conductivity is of order 2mS/cm , 
which is in the realistic range for neocortical pyramidal 
neuronssi. 

The input ISI of 4 — 8ms is important to understand- 
ing the dynamics of the Hodgkin-Huxley model. In the 
chaotic region between the k = 2 and A: = 3 locked states 
the coefficient of variation of ISI has a singularity at the 
transition between the odd-only and all-integer multi- 
ples of the driving period. The odd modes dominate 
in the vicinity of the /c = 3 state. The low-/c (high- 



frequency) bands vanish logarithmically near the line of 
critical points (gsyn^T^)- The firing rate has a minimum 
at T* + 0.2ms. Periodically stimulated giant axons 
of squid have similar nonmonotonic dependence of the 
firing rate on the current pulse amplitude between the 
k = 2 and k = 3 states^^. This experiment also showed 
linear dependence of the firing rate on pulse amplitude 
near the threshold for > T*, similarly to Fig. [TTJ 
Although the experimental pulses were rectangular, dif- 
ferent from the a{t) form with an exponential tail, the 
qualitative features do not depend much on the precise 
shape of a pulse. For short pulses the neuron's reaction is 
determined mainly by the time integral of the stimulus. 

The multimodal response occuring in certain sen- 
sory neurons may result from noise^^ or deterministic 
nonlinearity^^. It would be interesting to look for exper- 
imental evidence of the odd-all transition. It found, it 
would be a clear evidence that the neuron dynamics is 
dominated by nonlinear ity, not noise. 

The behavior of the model at small Ti may be useful to 
both coincidence detection and estimation of the signal 
strength. The optimal sensitivity in this case is inversely 
proportional to frequency. 

Our calculation also supports the view expressed by 
authors of Ref.^ that boundaries between various parts 
of the response diagram are not always clear-cut and may 
form complicated patterns. This statement also applies 
to the excitation threshold in the chaotic regime. 

In the Hodgkin's classification of intrinsic excitability^^ 
class 1 neurons maintain firing at arbitrarily low fre- 
quencies in response to weak inputs and have continu- 
ous frequency-current (/-/) curve. Class 2 neurons fire 
with certain relatively large frequency, usually of order 
40-50 Hz, when stimulus exceeds threshold and have a 
discontinuous /-/ curve. Class 1 and class 2 neurons 
sometimes are described as integrators and resonators 
respectively^^. According to the commonly held view 
a neuron cannot be an integrator and resonator at the 
same time. However we showed that the deterministic 
HH neuron in a chaotic regime near excitation thresh- 
old may oscillate with arbitrarily small frequencies and 
may perform integration at time scales much longer than 
the period of its main resonance. The character of the 
response depends strongly on the functional form of the 
stimulus and parameters of the model. A recent study 
showed that the same pyramidal neurons behave as inte- 
grators in vitro and resonators in vivo.^^ 

The multimodal response of the HH neuron near 140- 
180 Hz is not a typical resonance since no particular fre- 
quency is preferred. The multiples of the driving fre- 
quency alternate chaotically. The average output fre- 
quency depends nonmonotonically on the stimulus am- 
plitude. Similar nonmonotonic f vs. I relation was found 
in periodically stimulated giant axons of squid^^. Smaller 
stimuli favor higher multiples of the driving period. Stud- 
ies of large neuronal networks of various types suggest 
that there may be a complex interplay between the inte- 
grating behavior and the resonant action/^ 
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The ability to precisely control the nerve cell's po- 
tential oscillations is important in constructing de- 
vices performing the procedure known as Deep Brain 
Stimulation^-—, which operate at frequencies above 
100 Hz. While our model does not satisfy the charge- 
balancing constraint required in the stimulation of in- 
vivo systems, we believe the present study improves our 
understanding of high-frequency neural oscillators. 
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